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ABSTRACT 

We have constructed the mass function of globular star clusters in the Sombrero galaxy in bins of different 
internal half-mass density p/, and projected galactocentric distance R. This is based on the published measure- 
ments of the magnitudes and effective radii of the clusters by Spitler et al. (2006) in BVR images taken with 
the ACS on HST. We find that the peak of the mass function M p increases with p/, by a factor of about 4 but 
remains nearly constant with R. Our results are almost identical to those presented recently by McLaughlin & 
Fall (2007) for globular clusters in the Milky Way. The mass functions in both galaxies agree with a simple, 
approximate model in which the clusters form with a Schechter initial mass function and evolve subsequently 
by stellar escape driven by internal two-body relaxation. These findings therefore undermine recent claims that 
the present peak of the mass function of globular clusters must have been built into the initial conditions. 
Subject headings: celestial mechanics, stellar dynamics — galaxies: individual (M104, NGC 4594) — galaxies: 
kinematics and dynamics — galaxies: star clusters 
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1. INTRODUCTION 

One of the most fundamental properties of a system of star 
clusters is its mass function, defined here such that ip(M)dM 
is the number of clusters with masses between M and M+dM 
[i.e., ip(M) = dN/dM]. There is currently a debate about 
whether the shape of the mass function of old globular clusters 
primarily reflects dynamical evolution (Fall & Zhang 2001) 
or initial conditions (Vesperini et al. 2003). A key obser- 
vation here is that the turnover or peak of ip(M) at M p w 
(1-2) x 10 5 M Q is nearly "universal," varying little among 
galaxies of different masses and types and from one location 
to another within galaxies (Harris 2001). Interpreted naively, 
these facts tend to favor the scenario in which the present 
shape of ip(M) was imprinted when the clusters formed or 
were relatively young. On the other hand, the observed shape 
ip(M) w const below M p is a clear signature of long-term dy- 
namical evolution driven by internal two-body relaxation (Fall 
& Zhang 2001; Waters et al. 2006; Jordan et al. 2007). 

Another signature of relaxation-driven evolution, recently 
emphasized by McLaughlin & Fall (2007, hereafter MF07), 
is that the peak mass should be higher in subsamples of clus- 
ters with greater internal densities and thus higher evapora- 

1/2 

tion rates. The expected dependence is roughly M p oc p^ , 
where pi, = 3M/87rr| is the mean density within the median 
(i.e., half-mass) radius r/,. MF07 found essentially this depen- 
dence of M p on ph for the old globular clusters in the Milky 
Way, and they argued that this, together with the observa- 
tion ip(M) ks const for M < M p , provided compelling evidence 
for the scenario in which dynamical evolution determines the 
present shape of ip(M). 

In this Letter, we report on a similar study of the old globu- 
lar clusters in the Sombrero galaxy (M104, NGC 4594). Our 
goal is to check whether the MF07 findings apply only to the 
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clusters in the Milky Way or are representative of the clus- 
ters in other galaxies. We use the published measurements by 
Spitler et al. (2006) of the luminosities and effective radii of 
the globular clusters in images taken with the Advanced Cam- 
era for Surveys (ACS) on the Hubble Space Telescope (HST) 
as part of the Hubble Heritage Program (GO-97 14). The Som- 
brero galaxy is ideal for this project because it has a large ret- 
inue of clusters (N > 600), nearly all them are resolved in the 
ACS images, and there is relatively little contamination by the 
disk. 

2. OBSERVATIONS 

The ACS images are arranged in a 3 x 2 mosaic with 
overall angular dimensions of 600" x 400", corresponding to 
26 kpc x 17 kpc at the adopted 9.0 Mpc distance of the Som- 
brero galaxy. The images are moderately deep; for each point- 
ing, four integrations of 675 s, 500 s, and 350 s were made 
with the F435W (B), F555W (V), and F625W (R) filters, re- 
spectively. Spitler et al. (2006) selected candidate globular 
clusters based on the following criteria: ( i) integrated colors 
in the ranges 0.9 <B-R < 1.7 and 0.5 < B- V < 1.1; (it) ap- 
parent V magnitude brighter than 24.3 (roughly the 95% com- 
pleteness limit); (Hi) angular extent at least 0.3 pixels wider 
than the point-spread function (PSF); (iv) ellipticity less than 
0.5; (v) location outside the dust lane; and (vi) visual appear- 
ance more like a globular cluster than a star or galaxy. This 
produced a sample of 659 candidate globular clusters that is 
believed to be nearly complete brighter than My rj -5.5 and 
nearly free from contamination by foreground stars and back- 
ground galaxies. 

Spitler et al. (2006) present the magnitudes, effective (i.e., 
half-light) radii r e , and projected galactocentric distances R 
of all the clusters in their sample. 5 We estimate the mass of 
each cluster from M = Ly(M/Ly), with the total V-band lu- 
minosity Ly (after correction for Galactic interstellar extinc- 
tion) and the adopted mass-to-light ratio M/Ly = 1 .5 M Q /L@ 
(a typical value for old globular clusters in the Milky Way; 
see McLaughlin 2000). We estimate the internal half-mass 

5 The quantities denoted here by r e and R are denoted by Shi and Rqc in 
Table 2 of Spitler et al. (2006). 
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FIG . 1 . — Mass function of globular clusters in the Sombrero galaxy. The panels on the left-hand side show three ranges of internal half-mass density pi, in units 
of Mq pc~ 3 , as indicated (excluding clusters with R < 2. 1 kpc). The panels on the right-hand side show three ranges of projected galactocentric distance R in kpc, 
as indicated (excluding clusters with pi, < 30 Mq pc~ 3 ). The vertical bars represent Poisson (vW) errors. The solid curves were computed from equation (1) as 
described in the text by summing over terms with individual values of p;, , for all clusters in the pi, and R bins specified in each panel. The dashed curves, which 
are the same in all six panels, were computed from equation (1) with a single term with the median value of pi, for the full sample of 600 clusters. 



density from the definition ph = 3M/87rr| and the stan- 
dard 3D-2D conversion, r/, = (4/3)r e (neglecting any internal 
mass/light segregation; see Spitzer 1987). 

Figure 1 shows the main results of this paper: the mass 
function of the globular clusters in the Sombrero galaxy in 
three bins of internal half-mass density p h (left-hand panels) 
and three bins of projected galactocentric distance R (right- 
hand panels). We plot dN/dlogM rather than dN/dM for 
ease of comparison with the more familiar magnitude distri- 
butions. Based on the tests described below, we exclude the 
59 clusters with either p h < 30 M Q pc" 3 or R < 0.8' (2. 1 kpc), 
because the sample becomes less complete at low internal 
densities (for any R) and small galactocentric distances (for 
any ph). The mass functions plotted here show the familiar 
single-peaked shape. What is striking, however, is that the 
peak mass increases from M p w 1 x 10 5 Mq in the low- ph bin, 
to M p « 2 x 10 5 M Q in the middle-p/, bin, to M p « 4 x 10 5 M Q 
in the high-p/, bin. At the same time, the peak mass remains 
roughly constant at M p w2x 10 5 Mq for all three R bins. 
Both of these results have high statistical significance and are 
similar to the MF07 findings for the globular clusters in the 
Milky Way (see their Fig. 2). 

We have replotted Figure 1 with many different binnings in 
M, p^ and R, and always recover the same result — that M p 
increases with ph but is nearly constant with R — to within the 
statistical uncertainties. We have also made diagrams like Fig- 
ure 1 for subsamples of clusters divided by color at B-R = 1.3 
and find only minor differences. This is relevant because the 
globular clusters in the Sombrero galaxy, like those in many 
other galaxies, have a bimodal distribution of colors. It is be- 
lieved that the red, metal-rich clusters (B-R> 1.3) are some- 
what younger than the blue, metal-poor clusters (B-R< 1 .3). 
The spectra of several red and blue clusters indicate, however, 



that they are all old, with ages that differ by less than 3 Gyr 
from each other and from the ages of globular clusters in the 
Milky Way (Larsen et al. 2002). 

We have also made independent measurements of the mag- 
nitudes and effective radii of all the clusters in the Spitler et 
al. (2006) sample. Our procedure is similar to theirs, with 
only two minor differences: our measurements of V include 
r e -dependent aperture corrections and our measurements of 
r e are from the V-band images alone. Both sets of measure- 
ments are generally in good agreement, the mean offset (us 
minus Spitler et al.) and RMS scatter being only 0.06 and 
0.03 in V and 0.04 and 0.07 in logr e . These small devia- 
tions are at the levels expected for the slightly different pro- 
cedures. As another check, we have inserted artificial clusters 
with specified properties into the V-band images and mea- 
sured them in the same way as the real clusters. These ex- 
periments indicate that the uncertainties in V and r e increase, 
while the completeness of the sample decreases, for diffuse 
clusters (ph < 30 Mq pc~ 3 ) and those near the galactic center 
(R < 0.8')- We have imposed these restrictions on ph and R in 
most of our analysis (as noted above), but we have confirmed 
that the relations between M p , ph, and R are nearly identical 
in the full Spitler et al. sample. Finally, we have repeated the 
entire analysis using only our measurements of V and r e , and 
find results that are practically indistinguishable from those 
based on the Spitler et al. measurements. These checks give 
us confidence in our claim that M p has a significant depen- 
dence on ph but not R. 

3. MODEL 

In this section, we explore the implications of the results 
displayed in Figure 1 . As we have already noted, the peak of 
the mass function M p is a factor of about 4 higher in the bot- 
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FIG. 2. — Scatter plot of internal half-mass density pi, against projected 
galactocentric distance R for all clusters in the Spitler (2006) et al. sample. 
The vertical dashed line shows the maximum value of R accessible in the 
3x2 mosaic of ACS images. 

torn left-hand panel than in the top left-hand panel, while the 
median internal density ph within these subsamples increases 
by a factor of about 13. This is close to the expected scaling 

l li 

M p oc p h ' for relaxation-driven evolution, although the com- 
parison is necessarily crude because the range of ph included 
in each of the panels is large. The puzzle of why M p has lit- 
tle or no dependence on projected galactocentric distance R is 
explained by Figure 2, a plot of pi, against R for all clusters in 
the sample. Evidently, there is a huge scatter and only a mild 
(anti)correlation between p^ and R, thus accounting nicely for 
the weak dependence of M p on R despite the relatively strong 
dependence of M p on p h . (See Fig. 1 of MF07 for the analo- 
gous situation in the Milky Way.) 

To make a more precise and informative comparison, we 
follow the procedure developed by MF07, which in turn is 
based on the simple ("classical") model for the escape of 
stars from tidally limited clusters driven by internal relax- 
ation alone (see Section 3.2 of Spitzer 1987 and references 
therein). In this model, the mass of each cluster decreases 
linearly with time, M(t ) = Mo - p ew t, where Mo is the initial 
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mass and p ev = -dM/dt oc M/t± oc p h ' is the constant rate 
of mass loss (f r h being the half-mass relaxation time). For a 
population of clusters with the same internal density ph and 
age t, the present and initial mass functions are related by 
ip(M) = ip Q (M + p ev T) (Fall & Zhang 2001). This implies that 
the present mass function will have a peak at M p w p ev r, a 
low-mass shape ij){M) w const (forM < M p ), and a high-mass 
shape ip(M) « ipo(M) (forM > M p ). 

For the initial mass function, we adopt the Schechter (1976) 
distribution, ipo(M) oc M@ exp(- M/M c ), with a power-law ex- 
ponent 13 = -2 and an exponential cutoff at M c = 1 x 10 6 M Q . 
This choice is motivated by the observed power-law shape of 
the mass function of young clusters in the Antennae galax- 
ies, Vo(M) oc M~ 2 for M < 1 x 1O 6 M (Zhang & Fall 1999), 
and the fact that the mass functions of old globular clusters 



in many galaxies decrease more rapidly than a power law for 
M > 1 x 10 6 M Q (Burkert & Smith 2000; Jordan et al. 2007). 
We assume this mass function holds at some "initial" time, 
near the formation epoch but after a brief period of "infant 
mortality" caused by the activity of massive stars (photoion- 
ization, winds, jets, supernovae), since we are concerned here 
only with the long-term dynamical evolution of the clusters. 
In our model, this evolution flattens much of the power-law 
part of Vo(M) (for M < M p ) but preserves the exponential 
part. Thus, fortunately, most of our results are not sensitive 
to our assumed initial conditions. 

The mass function for a population of clusters with a range 
of internal densities is just the sum of the mass functions for 
the clusters of each density. From the relations above, we 
obtain a formula equivalent to equation (3) of MF07: 

V>(M) = A i r 2 exp[-(M+ A,)/M c ]. (1) 

i 

1 12 

Here, A; = {p, ty r)i oc p h j is the mass lost over the lifetime of a 
cluster, A, = A,/£ , 2(A,/M e ) is the required normalization con- 
stant for each term (£2 being an exponential integral), and the 
sum is over all clusters in the population. This model is based 
on the same dynamical evolution and initial conditions as the 
Fall -Zhang (2001) model, but the internal densities of clusters 
are no longer coupled to their orbits and hence to assumptions 
about the space and time dependence of the galactic potential. 
Thus, equation (1) is equally valid for spherical and triaxial, 
static and evolving potentials. 

We have computed dN/dlogM = (M/ 'log e)ip(M) from 
equation (1) and plotted the results in Figure 1. The coeffi- 
cient B in the mass-loss formula 

A; = B(p Ki /\tfM & pc- 3 ) 1/2 10 5 M Q (2) 

was determined as follows. We first fitted a model with a sin- 
gle term in equation (1) — an "evolved Schechter function" in 
the nomenclature of Jordan et al. (2007) — to the mass func- 
tion for the full sample of 600 clusters. The minimum x 2 oc- 
curs for A = 2.4 x 10 5 M Q and hence p ev = 1 .8 x 10" 5 M yr _1 
for r = 13 Gyr. Associating this rate of mass loss with 
the median internal density of clusters in the full sample, 
p h = 1.1 x 10 3 Mq pc" 3 , gives B = 2.3. The resulting "aver- 
age" mass function for the globular clusters in the Sombrero 
galaxy is plotted as the dashed curve in each of the six panels 
of Figure 1 . This is very similar to the corresponding func- 
tion derived by the same procedure by MF07 for the globular 
clusters in the Milky Way. 

We next computed the mass functions for the three ph bins 
and three R bins from equation (1) with a separate term for 
each cluster. In doing so, we used the same coefficient B = 2.3 
derived above for the whole sample to relate the individual 
values of A, and /?/, ,. The results are plotted as the solid 
curves in Figure 1 . These provide good (although not perfect) 
representations of the observed mass functions, with gener- 
ally similar shapes and peaks at about the right masses. In 
particular, the model M p increases by a factor of about 3 be- 
tween the low-p/, and high-p/, bins, while remaining almost 

constant for the different R bins. The first of these relations 

1 12 

is slightly weaker than M p oc p h ' because of the exponential 
cutoff at M c = 1 x 10 6 M Q . 6 

The observed mass function is slightly narrower or more 
sharply peaked than the model curves, especially in the left- 
hand panels of Figure 1 . One factor that might contribute to 

6 The simple scaling M p oc p~J 2 is exact in the limit M c 3> M p (MF07). 
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this is residual incompleteness for faint (i.e., low-mass) clus- 
ters, particularly in the lowest-density bin. Another is our 
simplifying assumption that all clusters have the same mass- 
to-light ratio. In reality, there will be a distribution of M/Ly, 
and this will broaden the mass function somewhat relative to 
the luminosity function [which we have adopted as a proxy for 
ip(M)], A more speculative possibility that would further im- 
prove the fits is that the cutoff M c in the initial mass function 
might increase slightly with pi, instead of remaining fixed as 
we have assumed here. In any case, it is clear from Figure 1 
that our simple, relaxation-driven model captures the essen- 
tial features of the observed mass function even without these 
adjustments. 

4. DISCUSSION 

Our findings for the globular clusters in the Sombrero 
galaxy corroborate those of MF07 for the globular clusters 
in the Milky Way. The Sombrero sample is larger (N = 600 
vs N = 146) but, because of the modest sensitivity to ex- 
tended sources and the restricted angular extent of the HST 
observations, it covers smaller ranges of internal density (6 < 
p h < 5 x 1O 4 M pc" 3 vs 3 x I0~ f < p h < 6 x 1O 4 M pc" 3 ) 
and galactocentric distance (2 < R < 15 kpc vs 0.6 < R < 
123 kpc). In both galaxies, the peak of the mass func- 
tion increases with internal density — approximately as M p oc 

1 12 

p h — but remains nearly constant with galactocentric dis- 
tance. These relations are consistent with each other because 
ph correlates only weakly with R. Thus, we now see that 
whether the peak mass appears "universal" depends crucially 
on one's perspective. In terms of galactocentric distance, it 
may be; but in terms of internal density, the more physically 
meaningful quantity, it certainly is not. 

We have interpreted these observations in the context of 
a simple, approximate model in which clusters form with 
a Schechter initial mass function and are subsequently dis- 
rupted by stellar escape driven by internal two-body relax- 
ation. Our model predicts present mass functions similar to 
the observed ones, including the correct dependence on both 
ph and R. We compute mass-loss rates and hence ijj(M) di- 
rectly from the observed p/,, and only then do we use the ob- 
served pi,- R distribution to reexpress ip(M) in terms of R. Our 
approach completely avoids having to specify the pi,-R distri- 
bution in the past and thus sidesteps the interesting but diffi- 
cult problem of explaining physically how the present ph-R 
distribution came about (which ultimately will require a de- 
tailed understanding of the formation and evolution of both 
galaxies and star clusters). Instead, we assume only that the 
present pf, of extant clusters are good guides to their past ph 

7 The strong radial mixing that occurs within evolving galaxies is readily 
apparent in movies generated from simulations of cosmological structure for- 
mation, although, to our knowledge, this has not been quantified in terms of 



and thus to their average mass-loss rates over a Hubble time. 
For a more complete discussion of these issues, we refer the 
interested reader to MF07. 

Most previous models link the internal densities and hence 
the mass-loss rates of clusters to their orbits in static, spherical 
galactic potentials (adopted for analytical simplicity). To ac- 
count for the observed weak radial gradient in the mass func- 
tions of globular clusters, these models are forced to have or- 
bital distributions with strong radial anisotropy at large galac- 
tocentric distances, possibly more than is allowed by obser- 
vations in the Milky Way (Fall & Zhang 2001) and certainly 
more than allowed in M87 (Vesperini et al. 2003). However, 
as Fall & Zhang (2001) have emphasized, these conclusions 
are primarily a consequence of the simplifying assumptions 
about the galactic potential, not the underlying premise about 
the disruption of clusters. In more realistic models, galaxies 
form and evolve by a hierarchical series of collisions, merg- 
ers, and other accretion events. Galactic potentials in this case 
are time-dependent and non-spherical, and the orbital energies 
and angular momenta of clusters are not conserved. The orbits 
and positions of clusters are rearranged many times and any 
initial gradients in M p or pi, are inevitably weakened if not 
eliminated, consistent with the observations plotted in Fig- 
ures 1 and 2. By comparison, the idealized models with static, 
spherical potentials overpredict gradients in M p and ph and 
therefore represent a limiting (i.e., "worst-case") scenario. 7 

Our results also have an important bearing on recent at- 
tempts to infer the initial form of the mass function of globular 
clusters. Vesperini et al. (2003) have advocated a gaussian- 
like or lognormal initial mass function with a built-in peak 
near 10 s M Q to account for the observed shape and weak ra- 
dial variation of the present mass function in the context of 
models with static, spherical galactic potentials (see also Ves- 
perini 2000, 2001). As noted above, this is not the only, or 
even the most natural solution of this problem. Our model 
provides a direct counterexample to these claims: an ini- 
tial mass function with a power-law shape (for M < 10 6 M Q ) 
evolves into the present mass function with a peak at M p w 
2 x 10 5 Mq and with the observed dependence of M p on ph 
and/?. 
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